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We consider the implications of the bond valence model (BVM) description of charge 
distribution in reduced ceria phases (Ce02-x)[l,2] to the models used to describe 
electronic and ionic conductivity in these phases. We conclude that the BVM is 
consistent with both the small polaron model (SPM) and the atomistic models 
which describe the electronic and ionic conductivities respectively. For intermediate 
phases, i.e., x ~ 0.3, we suggest the possibility of low temperature metallic conduc- 
tivity. This has not yet been experimentally observed. We contrast the BVM results 
and the conventional description of charge distribution in reduced ceria phases. 
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Introduction 

Ceria is a technologically important material with applications in high temper- 
ature electrochemical devices [31IHI51I61I7] . catalysis pi9]llOI)llj and oxygen gas 



sensors [T2]IT3] . It is the capacity of ceria to rapidly take up and release oxygen 

that makes it valuable for these applications. The process can be described by 
a reversible chemical reaction: 



The oxygen exchange controls both the electronic and ionic conductivities in 
reduced phases, Ce02-x, and both channels of charge transport are important 
in typical applications. This is illustrated by the example of a solid-oxide fuel 
cell where ceria forms the anode. The oxidation of the fuel takes place on the 
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anode and involves the abstraction of lattice O from the ceria. When an O 
atom is removed, a vacancy site results and two extra electrons are left in the 
crystal lattice of the anode. For the fuel cell to operate, two processes must 
ensue. First, the electrons move to the cathode through an external circuit, 
which relies on the electronic conductivity of the ceria. Second, oxide ions 
migrate from the bulk to the surface of the anode where the chemical reaction 
occurs, which requires a sufficiently large ionic conductivity. 

A description of the charge distribution in stable structures of Ce02-x may 
elucidate the microscopic mechanism of both the electronic and ionic conduc- 
tivities. For a long time, it was widely held that when an O vacancy forms 
in bulk ceria, the associated charge localizes on two of the nearest neighbour 
Ce sites [T^15|16II17|18] . We will call this the standard picture. Recent re- 
sults from atomistic simulations challenge this view, indicating that dopant 
ions prefer to occupy the second or third coordination shell of the O vacancy 
[T9|20f21|[22] . These reports are consistent with our own results from the BVM 
for slightly reduced ceria PPI2], which also contradict the standard picture. 

In Fig. [H we summarize the results of the standard picture and the BVM. In 
the standard picture, the two extra electrons localize in the nearest neighbour 
shell for all compositions of reduced ceria. The results of the BVM are different 
for different composition ranges of the reduced ceria phases (i.e. for different 
ranges of x in Ce02-x)- We roughly define 'low', 'intermediate' and 'high' O- 
vacancy concentration ranges according to a; < 0.2, 0.2 < a; < 0.3 and x ~ 0.5, 
respectively. In the low range, the BVM predicts that the two extra electrons 
localize in the next nearest neighbour shell; in the intermediate region, they 
delocalize in the nearest neighbour shell while in the high region, the result is 
the same as that for the standard picture. 

Except for fully reduced ceria, the standard picture and the BVM give quali- 
tatively different descriptions of charge distribution in Ce02-x phases. We aim 
to determine whether the difference in charge distributions results in a cor- 
responding difference in the transport properties of these materials. We thus 
study the electronic and ionic conductivities of the material using appropriate 
models of charge transport, and ask how the charge distribution will affect 
model predictions. 

In Section [l] we provide theoretical background on electronic and ionic con- 
ductivity in doped ceria. We apply these results to study experimental data on 
the electronic conductivity in the low and intermediate reduction regimes in 
Sections [2] and [3| respectively. The ionic conductivity is examined in Section 
m In Section in we conclude. 
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Fig. 1. A comparison between the standard picture and the bond valence model 
in describing the charge distribution in reduced ceria phases, Ce02-i:. In this 
schematic, we define 'low' as O vacancy concentrations up to just over x = 0.2, 
intermediate for x ~ 0.3 and high for x ~ 0.5. We compare the two descriptions 
with respect to where the excess charge is located and whether that charge is local- 
ized or delocalized. The main difference is that in the standard picture, the excess 
charge is localized in the first coordination shell for all compositions, Ce02-x- This 
picture agrees with the BVM analysis only in the high composition range. In the 
low composition range, the excess charge localizes in the second coordination shell. 
In the intermediate region, the BVM analysis implies a delocalized excess charge 
located in the first coordination shell. 



1 Theory 



Towards understanding general properties of electronic conductivity in re- 
duced ceria we start with a review of the electronic band structures of Ce02 
and Ce203, which correspond to extremal oxidation states of the Ce oxides. 
We then go on to describe the phenomenology of incoherent electron mobility 
in Ce02-a: phases and the relationship between carrier and electronic struc- 
ture. It turns out that the ionic mobility, described as a process of O vacancy 
migration, has similar properties to incoherent electron mobility. Thus, an im- 
portant question is how the two types of charge carriers might be distinguished 
experimentally. 
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Fig. 2. A schematic of the band structures of (a) Ce02 and (b) Ce203. Both Ce02 
and Ce203 are insulators. The band gaps and bandwidths are approximate and 
estimated from references j23ll24|25|26|27] for Ce02 and references |28|23|27] for 
Ce203. In each case, the bottom band, shaded black (blue in colour) is the valence 
band which is predominantly O 2p in character. The Ce 4/ states are shown in light 
grey (red). In the case of Ce02, these states are very narrow and centered ~ 6eV 
above the valence band. For Ce203, the 4/-manifold splits up into two sectors, a 
low energy and a high energy sector. The low energy sector is singly occupied and 
consists of two states per Ce atom while the twelve states in the high energy sector 
are all unoccupied. In both Ce02 and Ce203, the conduction band is shown in grey 
(green) and are mainly of Ce 5d character. When one considers these band structures 
of Ce02 and Ce203, it is not obvious that there could exist an intermediate phase 
such as Ce70i2, where low temperature metallic conductivity might occur. 



1.1 Electronic Band Structures of Ce02 and Ce20^ 



Ce02 and Ce203 are insulators for entirely different reasons. As can be seen 
from Fig. ga), Ce02 is a band insulator with a band gap of ~ 6eV. The 
valence band is of mainly O 2p character while the conduction band consists 
predominantly of the Ce bd states. Since there is a negligible probability for 
electrons to be thermally excited from the valence band to Ce bd states at 
relevant temperatures, there is no electronic conductivity in stoichiometric 
Ce02. In slightly reduced ceria, electronic conductivity involves the 4/ states. 
The latter exist within a narrow band and thus describe electrons that are ap- 
proximately localized on individual Ce sites. Electron hopping between these 
sites involves a phonon-assisted mechanism. This process is thermally acti- 
vated at temperatures below 1000 K or so (a value typical of working device 
conditions). 

In contrast, Ce203 is a Mott insulator [2SI23]- The Ce 4/ states are partially 
occupied with one electron per site. If the Coulomb repulsion between elec- 
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trons was taken to be weak, then band theory would predict the system to be a 
metaL However, the compact / states have a high on-site Coulomb repulsion, 
U ~ 6eV, which prohibits the double occupancy of a site. Since electronic 
conductivity would require that an electron hops onto a Ce site which is al- 
ready occupied, this material is an insulator with a band gap related to U. 
Fig. [2] shows that the band gap between the filled valence band and the con- 
duction band is ~ 5.5 eV and that between the partially occupied Ce 4/ states 
and the conduction band is ~ 3eV. Since U is larger than this energy gap, 
it is more likely for the electrons to be thermally excited into the conduction 
band than doubly occupying the local / orbitals. In any case, the number 
of free carriers is small at relevant temperatures. At low temperatures, the 
small polaron model (SPM) has been used to describe electronic conductivity 
in reduced ceria phases We now briefly review the main features of this 
model. 



1.2 The Phenomenology of Charge Carrier Mobility in Ce02-x Phases 

The electronic mobility can be equivalently described within two theoreti- 
cal formalisms: the SPM of Holstein [3CT]I31II29] and the Marcus-Hush theory 
[32 33J. To model the ionic mobility, we consider O vacancy migration or O 
self-diffusion in the Ce02-x phases |MII35l)36j . 

The fact that electronic mobility in reduced ceria is activated means that the 
mathematical description of electronic and ionic charge transport is similar 
(unlike in a metal, where there is no activation barrier to electronic mobil- 
ity). We take advantage of this correspondence below, while being careful to 
point out the significant differences between electronic and ionic transport 
that remain. 

We define a vacancy cluster consisting of the vacancy together with all the 
atoms up to and including the nnn Ce sites. Ionic conductivity involves the 
migration of the vacancy cluster on the O sublattice. The occupied 0-lattice 
sites (so called normal sites) and the vacancy sites are assumed to be ener- 
getically inequivalent with the difference in energy given by AiJ*^^, Fig. |3](a). 
Similarly, to describe electronic hopping on the cationic sublattice, we assume 
two types of Ce sites - the Ce^"*" and Ce*^^ which are energetically inequivalent. 
The energy difference between these sites is AH^^^, Fig. ^h). 

To mobilize a charge carrier requires that its associated energy barrier be 
overcome, so an energy of at least AH^^^, for electrons, or Aif*^^, for ions, 
must be provided by thermal fluctuations. In the case of ionic conductivity 
there is an additional barrier energy, AHl, associated with jumping to the 
saddle point at the boundary of the first coordination shell of Ce ions. Once the 
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Fig. 3. A schematic of the phenomenological description of the mobihty of charge 
carriers in Ce02-x phases. lUustration (a) shows the energies relevant to ionic con- 
ductivity. The shaded and empty circles represent an O ion on an O site and an O 
vacancy respectively. The two sites are both energetically and crystallographically 
inequivalent. The occupied site is lower in energy by Aff*^^. As the O ion moves 
to occupy the vacancy site, it experiences a potential shown by the dotted curve. 
There is a barrier potential given by AH^ which arises from the saddle point cre- 
ated by the layer of Ce atoms which the O ion has to cross to reach the O vacancy 
site. The solid curve shows how the potential changes as a result of the lattice 
relaxation, A//*, which follows the transfer of the O ion to the vacancy site. In (b), 
the case of electron mobility is shown. Here, the shaded circles are Ce^"*" sites well 
away from the O vacancy site. The empty circle is a Ce^'*' site which is assumed 
to belong to the vacancy cluster. All Ce sites belong to the point group. The 
potential experienced by an electron as it hops to the Ce^"*" sites is shown by the 
dotted curve. The Ce'^"'" site is lower in energy by AH^^g. The potential that results 
from the lattice relaxation after the electron has transferred to the middle Ce site is 
shown by the solid curve. This illustration shows that there is a trapping potential 
acting on the Ce^"*" site due to its proximity to the O vacancy. 

reaches the vacancy site, the lattice relaxes to the new configuration. The 
lattice relaxation has an associated energy, Aif*. Thus the mobility barrier of 
an O ion is given by the sum: Aif^ = AH^ + AH^. For electronic conductivity, 
no barrier energy analogous to AH^ exists so only the energy AH^^^ need be 
provided for the electron can tunnel to a neigbouring Ce site. However, as in 
the ionic case, once the electron transfer has occurred, the lattice relaxes with 
an energy given by AH^. Consequently, the mobility energy of an electron is 
given by: AH^ = AH^. 

For a mixed conductor with several charge carriers, the total conductivity is 
given by: 

a = ^ (Tj = ^ UjqjfXj (2) 

j j 

where aj, nj, qj and fij refer to the conductivity, density, charge and mobility 
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of the jth species of charge carrier. The O vacancies and/or electrons are 
bound to the O vacancy cluster so to obtain 'free' charge carriers, energy 
corresponding to at least Aif*^^ and Aif^^^ must be provided to the system 
for vacancies and electrons respectively. The thermally activated free carrier 
population is given by 

«, = «;exp(^] (3) 

where /cs is the Boltzmann constant and T is the temperature. For itinerant 
electrons, the mobility is given by: 

/^e = (4) 

m* 

where e and m* are the electron charge and mass and is the mean free 
time for electron scattering. For SPM and ionic conductivity, the mobility 
is diffusive and activated. It is related to the diffusivity, Dj, by the Nernst- 
Einstein relation [36] : 

u - ^ (5) 

A general expression for the diffusion coefficient, Dj, is obtained from random 
walk theory which assumes uncorrelated motion of the charge carriers. It is 
found that this expression has the same form for both SPM electronic [29] and 
ionic [Mfmii] conductivity given by: 

Dj = 7,(1 - c,)ay, exp (6) 

where the quantity, 7^, has a different meaning between electronic and ionic 
conductivity as will be described below. The quantity, Cj, is the concentration 
of normal sites, i.e., sites which are occupied by O atoms for ionic conductivity, 
and by Ce^"*" ions for SPM conductivity. Hence Cj = ^ with 9j and Nj being 
the numbers of occupied and total available sites of type j respectively. The 
attempt frequency, z/q is reported to be ~ lO^^s"^ for ionic conductivity [31] 
and ~ lO^^s"^ for SPM electronic conductivity [22] which corresponds to the 
optical mode vibration frequencies. 

Finally, the quantity, 7^-, for the case of ionic conductivity is given by: 

7. = g-exp(^— j (7) 

with z being the coordination number in the O sublattice of the site from which 
the O ion hops. For slightly reduced Ce02, z has the same value as it does in the 
unreduced phase so that z = 6. For the electronic case, the value of 7^ depends 
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on whether the electron transfer is adiabatic or non-adiabaticl30ll3H29l: 



: adiabatic electron transfer 

: non-adiabatic electron transfer 



Here, t is the electron hopping matrix matrix element and A is the reorgani- 
zation energy. 

For the motional free energy, we have already noted the energy terms that 
enter into the calculation of Aif^ for both ionic and electronic conductivity. 
From Eqs. [6| [5} [3] and Eq. |2] we get: 

^ = E^. = E(V^)exp(^) (9) 

where, 

Ei = AHi + AHl, (10) 

A- = 7,(g|/A:)c(l - c)a]nUe^^ |^ Ag^ + A^Ls ^ ^^^^ 

Small polaron (SP) electron transfer can be equivalently described in the 
Marcus-Hush theory which provides a general expression for the calculation 
of the rate of electron transfer 132]: 



~ h./A^^k^''''^ \ AXksT ) 

where k^t is the electron transfer rate or jump frequency (s~^), and AG^ is 
the Gibbs free energy of the electron transfer reaction. This model of electron 
transfer enables us to calculate t from some experimental data on mobility 
and activation energies. 

The electron transfer rate, ket, is related to the diffusion coefficient [38)13911401 : 



Dj = (13) 



Here, z is the number of equivalent nearest neighbour sites onto which the 
electron can hop and d is the jump distance. As we only have experimental 
data of mobilities for slightly reduced ceria, we consider the Ce sublattice to 
be that of Ce02 so that z = 12 and d = a\/2/2, where a is the lattice constant. 

The activation energy and the reorganization energy are related to each other 
by: 

..^'^ (14) 
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Fig. |4] is a simple schematic of how SP transport in Ce02-a; phases can be 
viewed in the Marcus- Hush picture (See also Fig. 7.8 in Mahan [H]). For 
AG° = 0, from Eqs. 12 13, we have the relationship required to calculate the 
hopping matrix element from the mobility: 



h(ksT)"{iyi, 



exp 



(15) 



1.3 Parameters for Conductivity in Ce02-x Phases 



We summarise in Table [T] both experimental and theoretical results for the 
main parameters in the conductivity equation, Eq. [9} in both the electronic 
and ionic cases. 



1.4 Signatures of Hopping Conductivity 



Table [T] shows that the motional enthalpy for electronic conductivity is > 
0.2 eV. The fact that the electronic conductivity is activated is consistent with 
incoherent (hopping) transport associated with localized charge. However, the 
activation energy is much less than the energy gap from the 4/ states to the 
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(a) (b) 

Fig. 4. A schematic for the small polaron model in Marcus- Hush theory, (a) Electron 
is hopping away from a site which is part of the vacancy cluster, (b) Electron is 
hopping between sites well away from O vacancy sites. The energy parameters are 
A - the reorganization energy, Ea - the activation energy and AG^ - the Gibbs free 
energy of the electron transfer process. The AG^ ~ for hopping between sites 
which are both well away from the O vacancy sites. 
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Table 1 

Some parameters for conductivity extracted from experimental and theory (Ref. 
|l2j) results in Ce02-a; phases. AH^ and AH^ are the motional enthalpies for 
electronic and ionic conductivity respectively, AH^^^ and AH^g^ are the association 
energies for the electron and O vacancy respectively and Ea is the sum of the 
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obtained for Sm-doped ceria: Ceo. ssSmo. 1501.925-5 (SDC15) 
2 obtained for Gd-doped ceria: Ceo.9Gdo.1O1.95_5 (GDCIO) 
^ obtained for Gd-doped ceria: Ceo.8Gdo.2O1.9_5 (GDC20) 

delocalized 5c? band in these materials which, e.g., in Ce02 is ~ 6 eV (compare 
Fig.|2|. 



SP transport can also be investigated by measurement of the Seebeck coeffi- 
cient of a material. In this case, the temperature dependence of the Seebeck 
coefficient is used to determine whether or not hopping mobility exists. We 
have noticed that mainly two models have been used in the literature to de- 
scribe the Seebeck coefficient of SP transport. The ffist one is more applica- 
ble to broad band semiconductors and shows a 1/T temperature dependence 
[5011^: 

S4'^){^^^) (16) 



where a is a coefficient of order unity. Eq. 16 is valid at low to moderate 
temperatures, i.e., < Es- At high temperatures, the Heikes formula, 

which is temperature-independent should be used [52|53] : 



S= - 1"^^ (17) 




Chaikin and Beni argue that Eq. 16 should not be used for narrow band 



materials as it does not describe localized states well 
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1.5 Decoupling the ionic from the electronic conductivity 



Both the Hall and Seebeck coefficients have the same sign as the charge carrier 
and so can be used to establish the type of charge carrier involved. A negative 
Seebeck coefficient was reported by TuUer and Norwick [29] implying that the 
electrons are the majority carriers for compositions where < x < 0.24 in 
Ce02-.. 

Another method to distinguish between electronic and ionic conductivities 
is impedance spectroscopy [H]. Here, a Nyquist plot of the the impedance 
spectrum gives the so-called 'half tear-drop' profile for a mixed conductor 
which arises from two charge carriers with different time constants being su- 
perimposed on the same curve. The results of Lai and Haile [IS] appear to 
corroborate the finding by Tuller and Norwick [29] that electrons are the ma- 
jority carriers in slightly reduced ceria. However, comparison of the two sets of 
results is not a simple matter as the former researchers used Sm-doped ceria 
whereas the latter used slightly reduced ceria. 



2 Electronic Conductivity in Slightly Reduced Ce02-x Phases 



Experimental work on electronic charge transport in the Ce02-x phases has 
been interpreted in terms of the small polaron model [^6p9] . Tuller and Nor- 



wick found a reasonable ffi to Eq. 17 for the Seebeck coefficient of the Ce02-x 
phases [2S] • Within the small polaron description the propagation of the charge 
carrier through the lattice is associated with a propagation of a lattice dis- 
tortion. There is no reason why the polaron must be confined to the vicinity 
of the oxygen vacancy. It was shown in Ref. [2] from the BVM analysis that 
most of the excess charge in the crystal of Ceii02o is localized on next near- 
est neighbour Ce sites from the O vacancy. In the same report, a prediction 
was also made that all the excess charge in CeeOn should be localized on the 
next nearest neighbour Ce sites. Thus, we see that the charge distributions 
obtained from the BVM do not contradict the SPM for small values of x and 
we can make a connection between the two pictures. The transport data can 



be used to obtain an estimate of the hopping matrix element, t, from Eq. 15 
using a = 5.411 A for the lattice constant we find: 

t = 0.006 eV (18) 

An estimate of t can also be obtained from the Harrison method (see Section [s]) 
which gives t ~ 0.1 eV. The bandwidth of the 4/ level in Ce02, W, calculated 
from density functional theory (DFT) is ~ 1.4 eV [23I2S]- The bandwidth is 
related to the hopping integral hj W = 16t so we have t ~ 0.09 eV in agree- 



11 



ment with the result from the Harrison method. In contrast, the value of t 
in Eq. 18 appears too low. However, since ne{T) is also thermally activated, 
the mobility cannot be accurately determined from fits over a small tempera- 
ture range. Moreover, correlation effects might be important in reducing the 
effective mobility. 



3 Electronic Conductivity in Ce02-x Phases at Intermediate Re- 
duction 



We now consider the electronic conductivity in phases of intermediate reduc- 
tion. We discuss the example of Ce70i2. Since the charge distribution obtained 
from the BVM suggests a delocalization of the charge on the Ce(2) sublattice 
[T|2] . we include calculations of the hopping matrix elements by the Harrison 
method [M]- The results enable us to discuss the question of charge delocal- 
ization in this crystal in a somewhat quantitative way. 

A key structural unit for discussing the electronic conductivity in this crystal 
is the divacancy cluster shown in Fig. |5| It results from the removal of two O 
atoms from the parent fluorite structure of Ce02. Two distinct Ce sites are 
obtained; the 5*6 site with two nearest neighbour O vacancies and the i site 
with only one O vacancy neighbour. The 5*6 site has been called the divacancy 
site \^5\ and it forms a shared corner between two coordination tetrahedra of 
the O vacancies as illustrated in Fig. [5j If we consider one of the O vacancies 
in the divacancy, the bond valence results indicate that the charge delocalizes 
among the Ce(2) sites that form an equilateral triangle. As different triangles 
are connected by relatively short Ce-0 bonds, it is plausible that charge is 
delocalized throughout the entire crystal on the Ce(2) sublattice. If this is the 
case, we can discuss the low temperature electronic conductivity of Ce70i2 
based on a simplified picture of a system of connected equilateral triangles 
each with a delocalized charge at one-third filling. We can disregard the Ce(l) 
sites because the electrons from the O vacancy do not occupy these sites. It is 
probably the case that a large activation energy is associated with the hopping 
of an electron onto Ce(l) sites. Fig.[6]shows a schematic of this simplified model 
where we have also included the direct Ce-Ce distances to give an indication 
of the relative separation between the equilateral triangles. We see that the 
distances between the triangles range from 3.66 to 4.09 A. The latter is the 
intra-triangle distance. Thus, inter-triangle distances are shorter or equal to 
the intra-triangle distances and thus the inter-triangle electron hopping matrix 
element may be comparable to the intra-triangle matrix element. 

The direct Ce-Ce distances are summarised in Table [2] along with the matrix 
elements for direct /-/ coupling between neighbouring Ce sites, tff. Table 
[2] shows that all but one of the matrix elements for inter-triangle electron 
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Ce(l) Ce{2) o(l| 0(2) OVacancy 

Fig. 5. Geometry of the Ce70i2 divacancy showing the symmetry at the corner 
shared by the two shghtly distorted tetrahedra. The divacancy cluster shown con- 
tains all the O atoms associated with the divacancy unit. Note that the Ce(2) sites 
are coordinated by 7 O atoms and therefore for each Ce(2) site, four Ce-0 bonds 
are not shown in this figure. The site symmetry of the Ce(l) site is Sq while that of 
the Ce(2) sites is triclinic. Estimates of the site valences from the BVM give +3.67 
and +3.21 for Ce(l) and Ce(2) respectively. 

hopping are larger than the intra-triangle value of 0.01 eV. The electron hop- 
ping matrix elements between lattice sites, tff were calculated from Harrison's 
method of universal parameters 



Since electrons may become delocalized on the Ce lattice by indirect hopping 
via oxygen sites we also should consider relevant Ce-0 bond distances. Fig. [7] 
shows the Ce-0 bond lengths associated with the bonding of the Cei site to 
its ten nearest neighbour Ce sites. The Ce-0 bonds are summarized in Table 
[2] together with the matrix elements for electron hopping between Ce sites via 
an O site, i.e., f-p-f hopping, tefj. Indirect hopping between two Ce(2) sites 
via an intervening O site involves a two-step process in which the electron first 
hops from the first Ce(2) to the O site and then from the O site to the second 
Ce(2) site. If t/p is the matrix element for hopping between a Ce(2) site and 
an O site, the overall matrix element for the electron hopping between two 
Ce(2) sites via an O site, t^ff, is then given by Eq. (19): 



teff = (19) 

where Ef — Ep is the energy gap between the Ce 4/- and the O 2p levels, which, 
in the calculation, we assume to be ~ 2 eV. The matrix elements for indirect 
electron hopping given in Table [2] are calculated from this formula [Mll56f57] . 
We note that all the matrix elements listed in Table [2] are for the a-a inter- 
action between the two orbitals as they are the most favourable. From Table 
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Fig. 6. Schematic of part of the Ce(2) sublattice of Ce70i2 showing the equilateral 
triangles connected to the Cei site. The boxed numbers are site labels which we 
reference in the notation CQgiteiabei- The numbers on the straight lines connecting 
the Ce sites are the direct Ce-Ce distances. The equilateral triangles are defined 
by joining together the three Ce(2) sites which are nearest neighbours to the same 
O vacancy. In this schematic, the reference triangle is highlighted in dark grey 
(purple) and the rest of the other triangles (grey (green)) are nearest neighbours 
to the Cei site on the reference triangle. Each equilateral triangle contains two 
electrons resulting from the formation of the O vacancy they coordinate. According 
to the bond valence results of I, the two electrons are delocalized among all three 
Ce(2) sites of the equilateral triangle. Thus, the three Ce(2) sites coordinating the 
same O vacancy can be described as an equilateral triangle at one-third filling. 

[2| we see that the te// matrix elements for intra-triangle hopping are equal 
and have the value O.lSeV. The same value is obtained for electron hopping 
from the Cei site to the Ce4 and Ce7 sites. An electron at Cei has an equal 
probability of hopping onto sites Ce2 and Ces (intra-triangle hopping) or onto 
sites Ce4 and Ce7 (inter-triangle hopping). The matrix elements for an elec- 
tron at Cei to hop to Ces are also similar. If these matrix elements are a fair 
indicator of the tendency towards charge delocalization, then they support our 
conjecture that the charge should be delocalized throughout the Ce(2) sub- 
lattice of Ce70i2. This delocalization of charge throughout the crystal would 
imply that Ce70i2 should exhibit metallic conductivity at low temperature. If 
this is true, then Ce70i2 should have the highest low-temperature electronic 
conductivity of all the reduced phases of ceria. 

We have not found any reports of low temperature electronic conductivity 
measurements of Ce70i2. But some results for Pr70i2 and Tb70i2 may be ap- 
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Fig. 7. Schematic showing the Ce-0 bonds in the Ce(2) sublattice. Except in the 
case of the highhghted triangle, only Ce-0 bonds connecting the Cei site to one of 
its ten nearest neighbour Ce sites. The highlighted triangle consists of the three Ce 
sites (including Cei) which are nearest neighbours to the same O vacancy. There are 
two electrons delocalized in this triangle so that it can be viewed as an equilateral 
triangle at one-third filling. The dotted Ce-0 bonds are not connected to the Cei 
site but have been included to show how all the three Ce sites in the equilateral 
triangle are connected by Ce-0 bonds. The white numbers on the Ce sites are site 
labels which correspond to those in Fig. [7] whereas the numbers shown on the lines 
connecting the various sites are the Ce-0 bond lengths in A. The bond lengths 
are summarised in Table |2} All the Ce sites in the Ce(2) sublattice are equivalent 
and each is seven-coordinated to O atoms and belongs to an equilateral triangle 
where the charge is delocalized as shown for Cei. The identification of the sites is 
as follows: black circles - Ce, grey (teal) - 0(1) and light grey (orange) - 0(2). 

plicable. In a study of the electrical conductivities of Pr02-x and Tb02-a;, Rao 
et al. [HH] found the highest conductivities for the compositions corresponding 
to Pr70i2 and Tb70i2. The authors also found that the Seebeck coefficient 
vanished at this composition. It was suggested that the seven-coordinated 
cations with an average charge of -1-3.5 would be responsible for the observed 
high conductivity. We, however, emphasize here that the results of Rao et 
al. are for electronic conductivity at temperatures ~ 200-600 °C whereas the 
bond valence results are for low temperature. Electronic conductivity at high 
temperature may be dominated by electron-phonon coupling effects which are 
absent at low temperature. 

We did not find any report on the electronic conductivity of reduced ceria 
at the exact composition corresponding to Ce70i2. The results of Tuller and 
Nowick [29] reproduced in Fig. |8]are the closest we got. 

The data in Fig. [8] which support the SPM shows an increase in the activation 
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Table 2 

A summary of the Ce-0 bond lengths from Fig. [7] and the direct Ce-Ce distances 
from Fig. |6j Matrix elements for electron hopping calculated according to the Har- 
rison method of universal parameters are included. The matrix element, t^ff, refers 
to electron hopping between the Cei site and its eleven nearest neighbour Ce sites 
via an O site whereas tff refers to a direct hopping by /-/ coupling 
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energy for the electron hopping with increasing x to about x = 0.25 which 
is the upper limit of the composition range measured. An extrapolation of 
these results to the composition of Ce70i2 {x = 0.29) would suggest a higher 
activation energy for electron hopping. However, it is not clear whether or not 
an extrapolation is a valid extension in the interpretation of these data. As the 
authors pointed out, it is expected that at this composition, an n—p transition 
will occur which may invalidate any extrapolative interpretation. Thus the 
situation regarding Ce70i2 appears somewhat inconclusive from the data on 
the SPM we were able to find. That low temperature metallic conductivity (i.e. 
band transport) should occur in Ce70i2 is unexpected from a consideration of 
the electronic structures of Ce02 and Ce203. These two phases are relatively 
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Fig. 8. Variation of activation energy from the high temperature region of the con- 
ductivity curve (1000 °C), with nonstoichiometry, x. Adapted from Tuher and Now- 
ick US]. 

well characterized and form the boundaries of the Ce02-x composition range. 
On this composition range, Ce70i2 is closer to Ce203 than it is to Ce02. Ce02 
is a band (charge transfer) insulator whereas Ce203 is a Mott insulator. The 
crossover from incoherent (hopping) to coherent (band) transport will occur 
at a temperature such that 5{T*) > e^a/h. 



4 Results from Atomistic Models for Ionic Conductivity in Ce02-; 
Phases 



There is a large body of literature on the study of ionic conductivity in doped 
ceria [55|HIS]|^61|62|63|64j . It is known that the local crystal environment of 
the O vacancies affects their mobility, so defect cluster models that simulate 
the local environment of the defect are employed. Below, we review some key 
results of the literature and consider whether they are consistent with the 
BVM. 

The mechanism of bulk ionic conductivity in doped or reduced ceria can be 
viewed either as oxygen self-diffusion or vacancy diffusion (of course, oxygen 
atoms move in the opposite direction to oxygen vacancies) under an applied 
electric field. The ionic conductivity is proportional to the product of the 
charge and concentration of the carriers and has an Arrhenius temperature 
dependence with an activation energy 



The process of doping Ce02 with a trivalent oxide, M2O3, can be described 
by: 

M2O3 ^ 2M'c, + 30o + Vo (20) 



where M'^g is a dopant ion on a Ce^^ site. Note that for M=Ce, Eq. 20 



IS 



equivalent to reduction of Ce02, Eq. [TJ For free O vacancies, the conductivity 



at low temperature is described by Eq. 10 with AiJ^^^ = AS"^^^ = 0, so the 



17 



ionic conductivity is a linear function of the O vacancy (and hence dopant) 
concentration provided that AS*^ and Ei are approximately constant. How- 
ever, experimental observations show that both Ea and a are highly nonlinear 
functions of the dopant concentrations, which means that a concentration de- 



pendence of El should be assumed if Eq. 10 is to be used 



To explain the relationship between E^ and dopant concentration, O vacancy 
cluster models were proposed. In these so-called defect associate models, the 
O vacancies are bound to the dopant ions and form stable associates of the 
the form (M'^^I/q)' and {{M'cJ^Voy. The latter is a neutral cluster which is 
only relevant at high dopant concentrations. Note that the superscript x is 
the standard label for this type of latter vacancy cluster , it should not be 
confused with the x in Eq. [l] (which never appears in superscript). In order 
for an O vacancy to contribute to the ionic conductivity, it must be ionized 
to give a free O vacancy: 

(M'^.V^-)' - M'^, + Vo or ((M'^J^ VoT ^ 2M'^e + V^o (21) 

The ionization of an O vacancy defect cluster costs an energy AHl^^, and 
entropy AS'^^^ so that the ionic conductivity in these clusters is described 



by Eq. 10 with non-zero enthalpy and entropy of association. AiY^^^ is the 
association or binding energy between the O vacancy and the dopant ion, 
M'fj^. Some values are given in Table [TJ Ionic conductivity measurements give 
information about Aif^^^ but no information about the types of defect clus- 
ters in the sample or their local atomic arrangement. To obtain this detailed 
information, various techniques have been adopted. Most results describing 
the binding energies of the O vacancies to the dopant sites have come from 
atomistic modelling [65P|68f 69170] . Our discussion here is limited to low tem- 
peratures (up to ~ 600 K) where the formation of defect clusters is almost 
complete [ifTT] . 

AHlgg can be calculated for a particular geometry of the defect associates and 
the results compared to the experimental E'^. An important early experimental 
finding was that the ionic conductivity of doped ceria depended strongly on the 
dopant ionic radius [72] . The local atomic arrangement in a defect cluster (and 
hence AHl^^) also depends on the dopant ionic radius [I2|7T] . Fig. |9] shows 
the trends in the binding energies calculated from atomistic modelling|68]. 
The most favourable arrangement in the (Ce'j^^V^') cluster has the Ce^"*" ions 
in the second coordination shell which is what we found for Ceii02o and what 
we expect for CeeOn. A similar study was recently done by Wei et al.jTOJ and 



the relevant part of their results is summarised in Fig. 10 



The results of Fig. 10 are for the composition Ce3oM2063 which for M=Ce^"'" 
is the supercell most often used in DFT studies [26|15fl6ll9fl8] . We have in- 
dicated in the figure where results of Ce203 are expected to lie based on the 
ionic radius of Ce^"*". The results of Wei et al. iTOJ can thus be understood as 
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1.60 



1st nearest neighbour 

Znd nearest neighbour 

3rd nearest neighbour 




0.65 0.75 0.85 0.95 1.05 1.15 

Cationic Radius (A) 

Fig. 9. Trends in the binding energies (from atomistic simulations) of an O va- 
cancy to a dopant ion located in one of the first, second or third coordination shell 
plotted as a function of the dopant ionic radius. The binding energy is calculated 
by subtracting the energy of the defect cluster from the sum of the energies of 
the individual components. Thus a positive binding energy indicates a favourable 
configuration for the defect cluster. Adapted from Minervini et al. [68]. 

corroborating those of Minervini et al. [68J. From atomistic simulations, Pryde 
et al. found that the binding energy for a vacancy cluster with two Ce^~^ ions 
in the first coordination shell was significantly less favourable than those with 
polarons in the second coordination shell (0.2 eV compared to ~ 0.4eV)[73]. 
They also found that the defect energetics of clusters involving In^"*" ions were 
qualitatively similar to those for the Ce^^ ions. The results of Pryde et al. 
on small polaron geometries are consistent with the BVM but contradict the 
standard picture. Similar results were obtained for doped Z1O2 by DFT meth- 
ods [H]. We, however, note that in contrast to these results, Deguchi et al. 
[75] concluded from their extended x-ray absorption fine structure (EXAFS) 
results that the dopant ions preferred to be in the first coordination shell with 
respect to the O vacancy. 



The higher binding energies for defect associates in which the dopant ions 
are in the second coordination shell have been explained in the following way 
[73 .68] : (i) the Coulomb interaction between the Vq and M'ce charged defects 
which favours the first coordination shell for the dopant ions, (ii) lattice relax- 
ation which primarily has to do with the relaxation of the Ce^"*" with respect 
to the Vq and M'ce charged defects around it. Because of its large positive 
charge, the Ce'^"'" ion prefers to relax away from the Vq site towards the M'^^^ 
site. This mode of relaxation is not possible if the M'^-^g ion is a nearest neigh- 
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-0.9 -I 1 1 1 1 1 

0.95 1.00 1.05 1.10 1.15 1.20 

Ionic radius (A) 

Fig. 10. The binding energies of various vacancy cluster configurations involving 
two dopant ions. The notation used for the defect cluster geometry is i-j-k where 
if the two dopants are designated as Ml and M2, then Ml and M2 are ith near- 
est neighbours (nn) and Ml and M2 are the jth and fcth nn to the O vacancy 
respectively (See insert in Figure.). The binding energy in this work is obtained by 
subtracting the total energy of the supercell with the maximum energy from that 
of a particular supercell. By this definition, a negative binding energy represents a 
favourable structure. We have indicated the position of Ce^"*" in the figure based on 
its ionic radius. The results suggest that based on this interpolation for the Ce^"*" 
ion, configurations 3_2_2 and 4_2_2 are the most favourable for arranging two Ce^"*" 
ions around an O vacancy. This is consistent with our BVM analysis of Ceii02o- 
Adapted from Wei et al. [70] . 

bour to the Vq site. Thus, this mode of lattice relaxation always favours the 
second nearest neighbour site for the dopant ion and (iii) a component of lat- 
tice relaxation driven by ion-size effects which favour the nearest neighbour 
position for small dopants. This is because the small ions prefer lower coor- 
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dination. For M=Ce^"'", effect (ii) is the most dominant and an energy gain 
of ~ 5 eV has been reported [US] which explains why the second coordination 
shell is preferred for these ions. Although we have not discussed the situation 
for high dopant concentrations where the neutral trimer ((M'(^g)2 V^')^ is ex- 
pected to dominate, Minervini et al have shown that the behaviour of these 
more complex systems is qualitatively similar to that of the (M'^^Vq) defect 
cluster 1681. 



We note that the atomistic models use the generalized Mott-Littleton method 
[76] which assumes an ionic description of the crystal lattice and the shell 
model for the ionic polarizabilities of the ions [77] . It is possible that the 
assumption of an ionic crystal for Ce02 might not be consistent with the 
results of bond valence calculations. Also, the models do not include vacancy- 
vacancy interactions which might be important at high vacancy concentrations 
as reported recently by Pietrucci et al. 



Additional data are needed for a comprehensive comparison with bond valence 
results. Ideally, one would also like curves of the binding energy as a function 
of the dopant concentration. The dopant of interest is Ce203 and since the re- 
duced phases of ceria considered in this paper are Ce70i2 (— > 3Ce02+2Ce203), 
Ceii02o 7Ce02 + 2Ce203) and CegOii 4Ce02+Ce203) the mole frac- 
tion of the dopant will need to cover the range 0.4 to 0.2. These are well within 
the capabilities of atomistic modelling. The range could be extended to lower 
concentrations using the 2x2x2 supercell typical of DFT calculations and 



also in Fig. 10 This classic case is interesting as it will provide an opportunity 
to compare the results from ab initio electronic structure calculations with 
those from atomistic modelling. 



5 Conclusion 



For slightly reduced ceria, the charge distribution obtained from the BVM 
is consistent with the description of the electronic conductivity by the SPM. 
For the Ce02-2: phases in the neighbourhood of a; = 0.3, we anticipate low 
temperature metallic conductivity or the highest low temperature electronic 
conductivity of all the reduced ceria phases. This conclusion relies on the es- 
timation of electron hopping amplitudes from Harrison's method of universal 
parameters. The detection of this metallic conductivity will discriminate be- 
tween the standard picture and the BVM desriptions of charge distribution in 
Ce02-a: phases. We have also considered the compatibility of the BVM charge 
distribution and the atomistic models used to describe ionic conductivity in 
these phases. We found that the atomistic models and the BVM give a consis- 
tent picture of the location of the charge. However, the data are quite limited 
and it is thus not possible to make a comprehensive assessment. 
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7 Graphic Abstract 



Is low temperature metallic conductivity possible in Ce70i2? We argue that 
since the distribution of the excess charge in this crystal can be described by 
a simple model of a network of equilateral triangles each at | filling, band 
conductivity may be expected. 
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